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Qv , We consider the modification of a single particle Schrodinger equation by 

the inclusion of an additional gravitational self-potential term which follows 

.•^ , from the prescription that the' mass-density'that enters this term is given by 

C^ . m\ip{f,t)\'^, where ilj{r, t) is the wavefunction and m is the mass of the particle. 

^^ ! This leads to a nonlinear equation, the ' Newton Schrodinger' equation, which 

O , 

- — , ! has been found to possess stationary self-bound solutions, whose energy can 

I , be determined exactly using an asymptotic method. We find that such a 

. . ! particle strongly violates superposition and becomes a black hole as its mass 



k> i approaches the Planck mass. 



C^ 



I. INTRODUCTION 

Inspite of tremendous success of quantum mechanics its interpretational aspects con- 
tinue to puzzle us, particularly in regard to measurements and extrapolation to macroscopic 
levels. The process of measurement involves a collapse of the wavefunction to one of the 
eigenstates of the operator being measured, and is not describable within the framework 
of the normal unitary evolution of the quantum mechanics. Further the quantum descrip- 
tion when extrapolated to macroscopic domain leads to paradoxical situations, which arise 
when we superpose macroscopically distinguishable quantum states. These two problems 



are related to each other, as the measurement process at some stage or the other involves a 
macroscopic apparatus and a result of the measurement by definition is a definite number 
resulting from a "pointer position". This precludes states which are superpositions of dis- 
tinguishable pointer positions. The pointer position is not predictable, and is found at one 
of the possible positions whose probability is dictated by its quantum wavefunction. The 
very nature of our classical description is thus incompatible with quantum superpositions 
and attendent probabilistic interpretation. 

There have been suggestions that gravity has a possible role in the fact that spatial 
superpositions of macroscopic objects are not seen. Our paper, which was stimulated by 
some recent suggestions by Penrose, is also an attempt in this direction. Central to this 
set of proposals is the idea that the wavefunctions that involve superpositions of spatially 
separated wavepackets should also have superpositions of the gravitational fields associated 
with the distinct mass distributions of the wavepackets. As has been lucidly explained by 
Penrose, there are basic difficulties in combining the covariance principle of general relativity 
and superposition principle of quantum mechanics. The difficulty is best seen in studying 
the superposition of states that are spatially apart. The gravitational field states involved 
in superposition require different space-times among which a point-wise identification is not 
possible. Under these conditions, it is impossible to define a unique time-translation operator 
and the very concept of stationary state. According to Penrose this leads to an uncertainity 
of energy, which makes such spatially superposed states inherently unstable. 

At another level, heuristically, if we assume that a particle may be localised at most 
in a size given by its Compton wavelength, then for a particle whose mass becomes of the 
order of the Planck mass or more and whose size is taken to be of the order of its Compton 
wavelength, the particle could be a black hole. This suggests that there is perhaps a problem 
with the QM of such particles. 

One way to model these effects may be to include a non linear mass dependent gravita- 
tional potential energy term in the standard Schroedinger equation, by postulating a 'mass 
density' given by m|-?/'(r, t)p, where ip{f,t) is the wavefunction and m is the mass of the 



particle. This term will automatically violate superposition. ( we cannot give any a priori 
reason for gravity to violate superposition and generate such a term in the Schroedinger 
equation). 

Such a nonlinear modified equation, called the Newton Schroedinger equation has a long 
history P|-p^ in the context of works that look at the infiuence of gravity on Superposition 
in QM. This equation has been used in several many body contexts, such as Hartree Fock 
calculations in plasmas as well as in astrophysics []T5|-|T7[1 . (It is also known as Choquard's 
equation [|l8l and the Schroedinger-Poisson equation [J19l . The existence of stationary 
solutions has been known (see for example [|18[) for sometime) In a preceding paper pO 



we looked at that the energy eigenstates of the equation and using an asymptotoic method 



could find the exact eigenvalue pOf.We review this in Sec. I, II and III. 

In Sec. IV (A and B) we take up the question of the spatial superposition of such 
localized stationary solutions and compute the energy difference between the superposed 
and stationary states , which may be identified with the breakdown of superposition. We 
show that this is related to an energy uncertainty calculated by Penrose using the the 
difference in the free fall accelarations corresponding to the the mass distributions of the 
two distinct spatial components of the superposition. This energy difference is used in analogy 
to Penrose's [l^ calculation to get a time for state reduction for an arbitrary mass particle. 

In Sec. IV (C) we first observe that our NS equation is valid only up to a limiting mass 
and not beyond. This follows from the fact that the bound state energy of the stationary 
state goes as m^ and thus overtakes the rest mass at a certain value of the mass causing our 
particle to be unstable. 

For masses just below the limiting mass, using our exponentially localised wave function, 
we can calculate the expectation value of the radius of our bound state particle of mass 
,m. On substituting this radius in the expression for the horizon parameter we find that 
the particle becomes a black hole somewhat below the limiting mass. We also find that 
superposition is strongly broken for such values of the mass. 



II. GRAVITATIONAL SELF-INTERACTION 

The implementation of the above set of semiclassical and nonrelativistic approximations 
amounts to saying that the particle experiences a self-gravitational potential, which arises 
from the gravitational potential energy due to the mass density given by m\ilj{f,t)\'^, where 
ip{f, t) and m are the wavefunction and mass of the particle respectively. Incorporating this 
interaction in the Schrodinger equation, one obtains the following equation. 

^j^diA^ = -—V^^(f, t) + \/(f)^(f, t) + mVG{r)i)ir, t) (1) 

at 2m 

where V{f) is the external potential acting on the particle and Vcir) is the gravitational 
self-potential arising due to mass density obtained from the wavefunction of the particle 
itself. Thus Voir) is given by, 

Vo« = -G/'!^^<%fd'n (2) 

J \r — ri\ 
and equivalently, 

V VG(r) = 47vGm\ilj{f, t) \^ (3) 

In order to get an idea of the magnitude of self-coupling, let us consider the external potential 
, V, to be the Coulomb potential, V = — e^/r. We can write the equation in a dimensionless 
form by introducing a free length a. In terms of this, the time is measured in units of 
r = 2ma^/h, and the energy is measured in units of e = h/r = fi^ /2me?. Let, 

t = rt; r = ar (4) 

E = eE; ij = ^/a^ (5) 

In these units, our equations take the form, 

'dt ^ -^^(^'^) ^^f "^("^^^^ + VG(rXf,t) (6) 



and 



V2i>G(r)=47rC|^(f,t)|2 (7) 

where C is a dimensionless coupling constant given by, 

C=^ (8) 

By choosing the length scale 'a'to set the coefficient of the Coulomb term to be unity, 
that is a = h? /2me^ we can immediately read off the value of the dimensionless parameter 

2 / 2 

, C = Gim? /e^ = —, where a is the fine structure constant. This clearly shows that the 

additional gravitational term we have added is relatively negligible compared to the coulomb 
interaction unless masses get rather large ; m of the order of m^i^fa . 

Before we examine the solutions of this nonlinear equation, it is useful to note the two 
conservation equations for particle number (or mass) and energy. These have been already 
noted in |^. The ffist equation that is easily established is for particle density, p{r,t) = 
ip*{r,t)ip{r,t). It obeys the usual continuity equation given by, 

| + V.4 = (9) 

In this and the following equations, we shall be working in units defined above, without 
explicitly putting tilde over dimensionless quantities. Next we consider the energy functional 

E(^(f,t)) , 

E = J[\Vi^{f,t)\' + V{f)mr,t)\' + ^VG{f)mr,m (10) 

The time-dependence oi E{ip{r,t)) satisfies, 

dE 



^^ + J d'rV.JE = (11) 

where Je is given in [^. Note that energy density is quartic in tp, and thus can not be 
expressed as an expectation value of an operator. We should also point out by calculating the 
expectation value of the momentum operator we find the consistency that our gravitational 
self potential does not contribute and therefore not give rise to a force or any spurious motion 

iJH- 



III. STATIONARY SOLUTIONS 



Solutions to this equation have been considered from a variational point of view in |^ 
and numerically in |2T] and more recently in [^. We now show that S-N equation admits 



stationary solutions , as found in [^| , of the form, 

^{f,t)=e~'^'(j){r) (12) 

The reason is that for solutions of this form the self-potential Vq becomes time-independent. 
0(r) obeys the equation, 



E<p(f) = -V^<p(f) + V(r)(P(r) -C f ^^^rfVi0(r) (13) 

J r — ri 



We now present an analysis of the equation in the absence of an external potential. 
The time independent equation also admits a variational interpretation, and can be 
regarded as an extremum of the functional, 

Hm] = - jdMvm? - - / /^^^P^rfv.rfv (14) 

with the normalisation constraint / |0(r)| d^r = 1. 

It is useful to record that with this constraint the energy functional above , under the 
scaling of r — ;> br, has a scaling behaviour that yields a minimum under the variation of b , 
indicating that the solution is stable with respect to scaling . 

We analyse one class of solutions of the free particle equation by doing an asymptotic 
analysis [0. This procedure allows us to obtain eigenvalues E for these solutions exactly. 



but the wave function is known only in the region of large r. The key observation is that 
if is taken to be an exponentially localised function about an arbitrary point, it gener- 
ates an asymptotic potential Vq which is a monopole potential C/r with corrections that 
decay exponentially with r, where r measure the distance from the point around which is 
localized. 

These remarks are best illustrated by taking 0(r) to be the ground state hydrogenic wave 
function (l/7rcr^)2e^'"/'^. The corresponding self-potential is. 



Voir) = -Cll - e--/-(i + i)] (15) 

Substituting the solution in Eq.(15) and keeping terms to order e~^^'^, one finds, 

^g-r/<r = (A _ i-)e-^/'^ - -e-'-/" (16) 

ar a^ r 

which yields value of E and a to be, 

a = ^ and E = -^ = -Cy4: (17) 

The asymptotic wave function and the potential gets exponential corrections that are 
exactly calculable while the eigenvalue is not changed ( see IQ for details ). (Note that 
here 'r' or a in units of a , the time in units of r = 2ins? /h and the energy in units of 

IV. PHYSICAL IMPLICATIONS 

In this section we consider some physical implications of the results obtained above. The 
above results show that the gravitational self-potential leads to a self bound state of linear 
extent of order ^. Again for a microscopic particle like the electron this means a localisation 
over enormous distances of order lO^'^A, with a binding energy of 10~*^ Ryd. Thus on scales 
much smaller than a, we do not expect the gravitational potential to affect things, and the 
usual quantum mechanics should apply. 

A. SUPERPOSITION 

In the context of the stationary states of the NS equation we can examine the existence 
of spatial superpositions or cat states in somewhat greater detail. Note that the centre of 
the stationary localised states discussed above can be chosen arbitrarily as the free particle 
equation is translationally invariant. More specifically, we calculate the energy of a cat-state 
consisting of two self-localised humps with separation much larger than a. The wavefunction 
of such a state is, 



0e(r) = [a0.(r)+/30,(|f-6|)] (18) 

There is a subtle difference between the locahsed states of our Newton Schroedinger equation 
and those of ,for example, an Hydrogen atom. In our case there is no separate centre of mass 
(CM) coordinate , whereas in the case of the H atom there is a localised stationary state for 
the relative coordinate and CM wavefunction is free to be in a plane wave or localised state. 

It would then be appropriate to identify our stationary wave function with a coordinate/ 
localised state of the CM of the Hatom. Since the H atom has very large mass compared to 
mass of the electron, localising the atom costs far less energy than localizing the electron. 
If we then neglect the former , we can effectively ( i.e. approximately) use the relative WF 
by itself to describe the localised atom. This we do. 

First consider the case of the usual Quantum Mechanics with linear operators which is 
based on the Superposition principle ( that is in the absence of the non linear gravitational 
interaction term). In this case if we identify (ps{r) as the usual lowest s wave state of the 
Hydrogen atom , that is a stationary state with energy eigenvalue Eq centred at the origin, 
then (psilr — b\) is the same state centred at b. If the separation \b\ » a, the scale of the 
wavefunction , the superposed state above has the same energy eigenvalue. 

Now consider the case in the presence of the non linear gravitational interaction term, 
and let (ps denote the stationary state solution derived above. When b >> a the function 
is normalized with a^ + /3^ = 1. Under this condition the energy of the superposed state is 
not the same any more. This is because i) the superposed state is no longer a solution of the 
non linear NS equation and ii) the gravitational potential energy between the lumps comes 
into play. In fact for asymptotic separation between the lumps it is always higher,given by: 

Ec = T-{a^ + (3^)Kn - 2a^/3^Ki2 (19) 

where T and Ku = Vq are the kinetic and potential energies of the single stationary state , 
and Ki2 is the gravitational potential energy between the humps. Since K12 is of order 1/6, it 
is much smaller than Ku (for asymptotic separations) and can be neglected. It is then easily 
seen that energy ofthe superposed state is higher than the energy of the single stationary 



state and in this asymptotic case it is clear that this originates not from the gravitational 
potential energy between the humps ,but from the normalization of the superposed state and 
the non linearity ( in its dependance on the wave function) of the gravitational interaction 

We specifically choose ol= P = 1/v^ to conform to the example considered by Penrose 
|I^. In this case for asymptotic separation we get a factor of 2 in the numerator , cor- 
responding to two lumps and a factor of 1/4 from the normalization , giving a potential 
energy equal to Kii/2 (neglecting the K12 term as we consider asymptotic separations ) for 
the superposed state, which is exactly half of that for the single lump stationary state. 

It is clear from the solution that the \E\ = 1^^/21 = |i^ii/2| ; this is also a simple 
consequence of the virial theorem. The energy difference, AV , between the superposed 
state and the stationary state is then given by , AV = \Kii/2\ = \E\. This is a measure 
of sorts of the breakdown of Superposition in the NS equation, as it is zero for the case of 
linear Quantum Mechanics. 

B. The Penrose conjecture on the time of state reduction of Superposed states 

Penrose [|l4l has argued that such superposed states are unstable due to energy uncer- 



tainty arising from the mismatch of the two space times due to the two lumps respectively 
that make up the superposition. 

More specifically a typical energy scale is obtained by Penrose by integrating the modulus 
squared of the difference of the free fall accelarations due to each of the two localised lumps 
that make the coordinate superposition of the stationary states. This energy turns out 
to be the gravtational energy of the difference of the mass distributions due to the two 
lumps. Penrose chooses to term this the gravitational self energy (GSE) of the difference. We 
will use the same acronym - GSE. 

A typical time of superposed state reduction is then obtained by dividing h by this energy 
in analogy with the decay of an unstable particle. These observations are general , for any 
given mass distributions -not specific to those considered above for the NS particle. 
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We shall now show that there is a close connection between the GSE computed by 
Penrose |T^ (and others) and the energy difference between the superposed state and the 
stationary state , AV , above. Assume that the ' mass density' for a particle governed by 
the NS equation is given by , Ps{r) = m|0s(r)p . 

On the other hand the normalized superposed wavefunction is 

0c(r) = -^[0^(r) + 0^(r)] (20) 

where \[p^~{^)Jm = (j)^,{r) =0^(1^1) 

and Jp'^g{f)/m = ^^^(r) = (f)s{\f — b\) 

The mass density for the superposed state is Pc{r) = m\(j)c{r)\'^ 

We can then write down the gravitational energy of the superposed state 

£,,Wr£^ (21) 

2 J J \f-ff\ ^ ' 

In the asymptotic limit , that is when |6| >> cr , there are terms in the expression above 
that are exponentially suppressed by the argument —\h\r and can be dropped. In the neglect 
of these terms and with some routine juggling we can get the energy difference between the 
superposed state (pcij) and the stationary state (psij), AV = Vc — Vg 

^V=^-^l I ((^Hr))(p^(rO) + {p\r)){p\r')) - {p\r)){p\r')) - {p\r)){p\r')) ^,^^,^, 
42 J J If-r'l 

(22) 



which is proportional to the GSE A that appears in [|14| 



^ = ^^G f f ((/^'("))(/^'("0) + (P'(^))(P'(^0) - ip'ir))ip\r')) - {p\r)){p\r')) ^,^^,^, 
J J \r — -PI 

(23) 

This establishes the connection between the energy difference, AV^ ( which is the the 
difference of the gravitational energies ( GSE ) between the the superposed (asymptotic) 
and stationary state ) that is a measure of the breakdown of the superposition principle for 
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us and the energy uncertainty, A, (which is propotional to the GSE of the density differences 
of the two humps ) associated with superposition in ||14|| . 

There are problems with the above identification. We should note here that it is ob- 
vious that the mass density expression used in constructing the gravtational term for the 
NS equation is not in any sense mean that the electron mass is distributed according to 
its wavefunction it can have only the usual QM probalilistic interpretation. However as the 
mass of the particle is increased its behaviour becomes more and more classical and such 
an interpretation may become plausible. ( Anyhow, we expect that if the asymptotic sepa- 
ration parameter b is of the order of laboratory scales of centimeters ,and we constrain the 
stationary wavefunction to be spatially smaller ,say 10~^cm. This implies that realistically 
we are dealing with particles of a mass in excess of 10^'^Gev.) 

For a particle of mass ,m , governed by the NS equation we importantly find that ,AV = 
\E\, goes as m^ 

This is indeed different from the classical expression for the GSE of a constant density 
mass distribution which goes as ,m?/R , which in turn goes as m^^^ , for composite objects 
like droplets considered in [^ 



(We observe that these expressions for , AV^ and A are proportional to the GSE for the 
single statioary state because, for asymptotic separation, the gravitational potential energy 
between the lumps can be neglected. ) 

On assuming Penrose's conjecture that this is the energy uncertainty associated with an 
unstable particle, we can now explicitly find the time of state reduction for the superposed 
state of a NS particle as a function of m . This is 

Tsp = h/AV = h/\E\= 2h^/G'^m^ 

For specificity this is 

i)10™ sec, for an electron 

ii) 10~^° sec. for a particle of a mass that is less than the limiting mass by a factor of 5 
-which puts it it safely in a non relativistic , stable and non black hole regime ( see following 
section ). 

11 



For comparision we can estimate the time of state reduction for a droplet of constant 
density ,when the GSE given in |2^ goes as , rrv'/^. 

Tsp = h/A = h/Gmyn = h/Gm!'/^ 

This is 10~^° sec. for a particle of a mass that is less than the limiting mass by a factor 
of 5. 

C. LIMITING MASS OF A PARTICLE 

It must be noted that the energy eigenvalue ( which is negative for the bound state ) goes 
as TJV' . This results in an instability as the sum of the rest mass energy and the attractive 
binding energy becomes negative. At this point our Shroedinger Newton description with 
only the gravitational interaction breaks down. The value of this limiting mass is obtained 
on putting the ratio of the modulus of the energy eigenvalue and the rest mass , F = 1 (Note 
F must be small compared to unity for our analysis to work.) 



,2h\^ 



1 



m„ = (-^)4 (24) 

This is effectively the Planck mass (apart from a constant factor of 0(1) ) Perhaps this is 
not so surprising for a theory which has only gravitational interactions. 

We would like to consider more of the physics as we approach this limiting mass. To 
this end we determine some relevant parameters for particles of mass that approaches the 
limiting mass and compare them to their counterparts for the electron. We note that our 
theory has only one dimensionless parameter, 5 = {m^G/hc). All dimensiomless quantities 
are then expressed in terms of, 6. 

i) The parameter F = \E\/mc^ above goes as ^irn^G/TicY = ^6^. It is 

1 for the limiting mass , 

10"^*^ for the electron 

Recall that in our theory the breakdown of Superposition is linked to the energy , AV 
which is further equal to , E , the energy eigenvalue of the bound state. Dividing by the rest 
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mass energy yields, F = \E\jm(? , which then is also a measure of sorts of the breakdown 
of superposition . 

ii) The size parameter R 

The typical size for a particle, R ,may be evaluated by calculating the expectatation 
value of , 'r' , for the wave function (l/Q;^7r)2e~'"/",where, 

(x = ao = a2/C = Ti^ /Guv" 

Thus , i? = f 

iii) The Horizon parameter , 2Gm/ Rc^ = ^{Gm^ /TicY =|5^ 

where we have substituted for R above . We should further notice that this parameter 
goes as m^ or the square of 5. 

This parameter decides if a mass distribution is a black hole or not. Specifically, if it is 
greater than 1 we have a black hole. It is 

I for the limiting mass 

10^^° for an electron 

Clearly, at the limiting mass the particle is a black hole ! 

This analysis shows that in the presence of a gravitational modification of the 
Schroedinger equation as given above we find not only a stability problem with masses 
as we approach the limiting mass but that such a particle would be a black hole. Note 
however that a reduction of the mass by only a factor of 5 of the limiting mass gives a very 
sensible non relativistic description and brings down the Horizon parameter by a factor of 
125 taking it safely away from a black hole. 

We have studied a certain cojectural non linear modification of the non relativistic 
Schrodinger equation due to gravity - the Newton- Schrodinger equation to look at the 
quantum mechanics of such a particle. The implication is that in this description as we 
approach the Planck mass not only do we strongly violate superposition but the particle 
becomes a black hole, thereby putting a limit on the mass of elementary particles. 
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